Temperature suppression of STM-induced desorption of hydrogen on Si(lOO) surfaces 
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The temperature dependence of hydrogen (H) desorption from Si(lOO) H-terminated surfaces by a 
scanning tunneling microscope (STM) is reported for negative sample bias. It is found that the STM 
induced H desorption rate {R) decreases several orders of magnitude when the substrate temperature 
is increased from 300 K to 610 K. This is most noticeable at a bias voltage of —7 V where R decreases 
by a factor of 200 for a temperature change of 80 K, whilst it only decreases by a factor of 3 at 
—5 V upon the same temperature change. The experimental data can be explained by desorption 
due to vibrational heating by inelastic scattering via a hole resonance. This theory predicts a weak 
suppression of desorption with increasing temperature due to a decreasing vibrational lifetime, and 
a strong bias dependent suppression due to a temperature dependent lifetime of the hole resonance. 
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Desorption of hydrogen (H) from a H-tcrniinatcd sili- 
con (Si) surface has recently been studied actively [|l|-[lO[ , 
because of its potential applications to construct novel 
optical and electronic devices. For positive sample bias, 
Shen et al. g have shown that below 4 V the STM 
induced H desorption is related to vibrational heating 
0,0 by inelastic scattering of tunneling electrons with 
the Si-H 6ct* resonance, and above 4 V it is related to 
direct electronic excitation of the Si-H bond. At nega- 
tive sample bias, the desorption mechanism has recently 
been explained by the authors |1^ to be due to inelastic 
scattering of tunneling holes with the Si-H Scr hole reso- 
nance and quantitative agreement between experimental 
desorption data and first principles calculations was ob- 
tained. 

The present paper focuses on the temperature(T) de- 
pendence of the STM induced H desorption rate(i?) from 
Si(100)-H(2xl) surfaces in ultra high vacuum (UHV) at 
negative sample biases from —10 V to —4 V and temper- 
atures from 300 K to 610 K. For a given sample bias(Vb) 
and tunnel current (/) we measure a decreasing R with 
increasing T. The effect is most pronounced at Vb = — 7 
V, where R decreases by a factor 200 for an increase in 
T of 80 K, whilst it only decreases by a factor 3 at Vt, 
= —5 V upon the same temperature change. Previous 
studies of the T dependence of R have been carried out 
for positive sample bias and at T = 11 K, where R was 
reported to increase by a factor of ~ 300 compared to 
room temperature |9). The lower R with increasing T 
was explained by a decreasing vibrational life time of the 
Si-H stretch mode with increasing T ||l^,|j| . 

In the present study the T dependence of R has been 
modelled using a first principles model |15|| of inelastic 
scattering of tunneling holes with the Si-H Scr hole reso- 
nance. Two temperature effects, which are based on the 
T dependence of the lifetime of vibrational excitations 
and electronic excitations, respectively, are included in 



the model. At Vb = — 5 V, the T dependence of R is re- 
lated to the lifetime of the vibrational excitation, which is 
similar to the result for positive Vb reported previously 
by Foley et. «/.§]. At 14 = -7 V the T dependence 
from this effect is too weak to explain the experimental 
data. The strong T dependence of R in this case is mainly 
caused by the T dependence of the lifetime of the ba hole 
resonance, which gives rise to an additional strongly bias 
dependent temperature effect. The combination of a T 
dependent vibrational lifetime and a T dependent elec- 
tronic lifetime accounts very well for the experimental 
data in the whole bias range from —10 V to —4 V. 

The experiments were performed at a base pressure of 
- 7 X 10"^ Pa using n-type Si(lOO) {No = 1 x lO^^ cm^^) 
samples and a JEOL JSTM-4000XV microscope with 
electrolytically sharpened tungsten (W) tips. Samples 
were heated using a direct current flowing between two 
contacts at each end of the samples, the STM tip being 
positioned in the middle of the two contacts. A change 
in Vb caused by the current flow was corrected for in 
the STM bias circuit, and the absolute value of Vb was 
verified by measuring the zero point in current-voltage 
curves. The Si(100)-H(2xl) surfaces were prepared by 
standard procedures [0 . The surfaces are constructed by 
an array of parallel 'dimer rows', where the dimer con- 
sists of two dimerized Si atoms each bonded to a single 
H atom. If one of the H atoms is desorbed, a Si dan- 
gling bond appears and is observed as a bright spot in 
filled-state STM images. The desorption rate (i?) of H 
atoms is determined by the following STM parameters: 
the scanning speed (s) of the STM tip, Vb, / and T. The 
experiments were carried out at negative Vb between —10 
V and -4 V and at T between 300 K and 610 K. Increas- 
ing T above 610 K caused an excessively large thermal 
desorption of H atoms and the study is therefore limited 
to T < 610 K. 

In order to establish the relationship between i?, /, 



Vb and T, the surface was scanned by the STM tip at 
various conditions, and R could be determined from s 
and the number of desorbed H atoms identified from the 
STM images. Two different kinds of experimental meth- 
ods were employed. In one method the STM tip raster 
scanned areas of 20 nm x 20 nm on the Si(100)-H(2xl) 
surface using a fixed value of Vb, while / and s were ad- 
justed to yield a suitable countable number of desorbed 
H-atoms (typically 100-200 atoms). This method was 
used to determine R as function of / for constant Vb- In 
another method, line scans with a fixed s = 2 nm/s were 
performed at various Vb, and / was adjusted to yield des- 
orption of 50% of the H atoms along a line scan on the 
surface. This method was used to determine the rela- 
tionship between Vb and / for constant R. 

Figure 1 shows a filled state STM image recorded at 
T = 530 K after desorption of H from the Si(100)-H(2x 1) 
surface. The bright line running from the top right to the 
bottom left is a chain of Si dangling bonds created by 
scanning the tip along a line with s — 2 nm/s, Vb — —7.0 
V and 1 = 8.0 nA . The width of the line is approximately 
equal to the width of one dimer row and the resolution 
is similar to that obtained for lithography at room tem- 
perature po[ |. 

Figure 2 shows R plotted as function of / for (a) 
Vb = — 7 V and (b) Vb = —5 V and for temperatures 
T = 300 K (squares), 380 K (crosses) and 450 K (tri- 
angles). We have made least squares fits to the data of 
a power-law dependence of R upon I, R = i?o(-^/^dcs)", 
where /des is the tunnel current corresponding to i?o = 4 



At Vb 



-7 V the least-squares fitted parame- 



ters are: a = 5.7 ± 0.7, 4.7 ± 0.5 and 4.6 ± 0.7, and 
/dos = 0.85 nA, 2.0 nA, and 3.3 nA for T=300 K, 380 
K, and 450 K, respectively. At Vb = —5 V, the corre- 
sponding parameters are: a — 6.3 ± 1.3, 4.6 ± 0.5 and 
6.4 ± 0.9, and /dcs = 4.5 nA, 4.2 nA and 6.4 nA. The 
data follow the power law very well at all temperatures 
and the slope is almost independent of T. The tempera- 
ture effect is largest at —7 V; for example at / = 2 nA, 
R decreases from 700 s^^ to 3 s^^ when T is increased 
from 300 K to 380 K, while at -5 V and / = 8 nA, i? 
only decreases from 120 s^^ to 40 s^^ upon the same 
temperature change. 

The relationship between /dcs and Vb, corresponding 
to a fixed desorption rate of i?o = 4 s^-'^ or 50% desorp- 
tion of H along a line scan with s = 2 nm/s is depicted 
in Fig. 3 for T ranging from 300 K to 610 K. At a given 
T and Vb, i? > 4 s^^ for values of / > /dcs- At all T, 
the data show a minimum in Ides at Vb '^ — 7 V, and 
therefore a maximal desorption yield at this voltage. 

Previously we have shown that the room temperature 
data can be explained by vibrational heating of the Si-H 
stretch mode due to inelastic scattering with the Si-H 5cr 
hole resonance, and first principles calculations of R due 
to this mechanism are in quantitative agreement with 
the experimental data |1^. In the following we use this 
method to calculate the T dependence of R and compare 
it with the experimental data. 



The first principles method for calculating the inelastic 
current is based on a high voltage extension [|6[ of the 
Tersoff-Hamann model 11^ for the STM tunnel junction, 
and has been described in detail elsewhere |l^,|l^. We 
include inelastic scattering events with energy transfer 
ibTiwoj ±2?iwo and ±3?iaJo: where hoj^ = 0.26 eV is the 
energy of the Si-H stretch mode. Energy relaxation due 
to coupling with silicon phonon modes is described by 
a current independent energy relaxation rate, 7t, which 
we determine from experiment pq] in a manner similar to 
Ref. P] . To obtain R, we solve the Pauli master equation 
for the transitions among the vibrational levels of the H 
potential well, and assume that desorption occurs when 
the energy of the H atom exceeds the desorption energy 
Edcs = 3.36 eV, corresponding to a truncated harmonic 
potential well with 13 levels. 

In the calculation of the inelastic current a correction 
is made for the difference in the Fermi level, £|, and the 
band bending, $, between the slab and the sample. This 
is carried out by using an offset between the sample bias 
of the slab model Vb and the experimental Vb 



eVb = eVb+eJ,-<l-(£^-$), 



(1) 



where the values with a tilde are slab quantities. Both 
e'^p and $ depend on the sample temperature. However, 
since the changes in $ and e|^ nearly cancel each other, 
the offset, Vh — Vh — 0.30 ± 0.05 V, is almost constant 
for the temperatures considered in this work. The de- 
pendence of the inelastic current on T due to changes in 
$ and e^p can therefore be neglected. 

We first consider the effect of a T dependence of 77- . 
The relaxation of the Si-H stretch mode due to phonon- 
phonon interaction is expected to involve three quanta 
of the Si-H bending mode, hub — 0.078 eV, and a Si 
phonon n^si = 0.026 eV ^^. The T dependence of 
this interaction is described by fl3] 

It = 7o{[l + nT{ijJh)f[^ + ?^T(^Si)] - nT(t^b)^"-T(wsi)}, 

(2) 

where nriuj) = l/(exp(?icj/fcT) — 1) is the Bose-Einstein 
occupation number. Using the experimental room tem- 
perature relaxation rate, 7300 = 10^ s""'^ and Eq. (0) 
the relaxation rate at T = K is determined to be, 
7o = 5 X 10^ s~^. At 600 K the relaxation rate obtained 
from Eq. (||) is , 7600 = 3 x 10^ s~^, which is 3 times 
larger than the value at room temperature. 

We note that lateral diffusion of the excitation into the 
H overlayer through incoherent exciton motion is not im- 
portant for the present system. This is due to the fact 
that inelastic scattering with an energy transfer of 2hujQ 
dominates the desorption process and the lateral diffu- 
sion rate for excitations with n > 1 is very low due to 
the anharmonicity of the Si-H bond potential |jlj,p3 ■ 

The dashed curves in Fig. 2 show the result of this 
model for the dependence of i? on / for Vb = —7 V 
and Vb = -5 V and T = 300 K, 380 K and 450 K. 



We note that there are no adjustable parameters in the 
model. For the present system inelastic scattering events 
with energy transfer 2fiuJo give the dominant contribu- 
tion to R in the temperature and voltage range investi- 
gated. Since there are 13 levels in the truncated harmonic 
potential well of the H atom [gO| it can be shown that 
R ~ (h/^/TY^^'^ 24 1 where I2 is the inelastic current with 
energy transfer 2hLdQ. From Eq. (^ we have 7t = 1, 1.3, 
and 1.7x lO^s-i at T = 300, 380, and 450 K, respectively. 
From this model we obtain iT!(300)/i?(380) - 1.3^-5 ~ 5 
and i?(300)/i?(450) - 30. At 14 = -5 V this is in good 
agreement with the measured T dependence, whilst the 
measured T dependence is almost two orders of magni- 
tude larger at Vb = — 7 V. The dashed curves in Fig. 3 
show the calculated bias dependence of /dcs for temper- 
atures, 300 K, 380 K, 450 K, 570 K, and 610K. The 
minimum in /dcs occurring at —7 V is reproduced by the 
model and is related to the position of the energy level of 
the 5ct hole resonance iQ. For all voltages R decreases 
with T. For Vb > — 5 V there is quantitative agreement 
between the theoretical model and the experimental data. 
However, for Vb < —5 V there is an additional tempera- 
ture dependence not accounted for by the T dependence 
of 7t. 

As already discussed for the present system there are 
no significant temperature effects due to band bending 
and lateral diffusion of the vibrational excitation, and 
calculations of the tunnel current using the appropri- 
ate Fermi-Dirac occupation factors show that the ther- 
mal current is negligible. At the elevated temperatures 
we still have atomic resolution for imaging and lithog- 
raphy(see Fig. 1), implying that the temperature does 
not affect the tip shape. At elevated temperatures the 
deposition of H atoms from the tip to the surface be- 
comes important pll , but this effect is not large enough 
to change i? by a factor of 100. 

Instead we consider the effect of a T dependent life- 
time of the 5(7 hole resonance. The time evolution of the 
wave function, (p^a, is given by 






ient/?ig-t/2Tph 



(3) 



where ^/i^ are eigenstates of the adsorbate-t-substrate sys- 
tem with eigenvalue, e^, and due to electron-phonon cou- 
plings they have a finite lifetime, Tph, which we take to 
be independent of /i. If the 5cr resonance is approximated 
by a Lorentzian line shape we obtain 



05a(t) = 05a(O)e- 



-i£5^t/''p-(l/2-ro + l/2rph)t 



(4) 



where, £50-, is the resonance energy and To is the life- 
time due to electronic coupling with substrate eigen- 
states. The lifetime of the resonance may now change 
with temperature due a T dependence of both Tph and 
Tc. The T dependence of Tph for different eigenstates of 
Si has been measured in Ref. [g2|,|3|, and ?i/2Tph '--^ 0.1 
eV at room temperature and h/2Tph ~ 0.2 eV at 600 K 
[Ea,E3|. We also expect a T dependence of To, since the 



electronic coupling between H and the Si surface depends 
on the position of the Si atoms. 

The energy dependent probability of inelastic scat- 
tering, Pn{e), with transfer of energy, nhcu, is for a 
Lorentzian resonance given by p4[ 



P„(£,At)cx 



A. 



[{e 



£50 



A2,]"+l^ 



(5) 



where At = l/2Tph + 1/2tc is the Half Width at Half 
Maximum(HWHM) of the resonance at temperature T. 
From P„ the T dependent inelastic currents can be cal- 
culated using 



/„(e,T)«/„(e,0) 



P»(£,At) 
Pn(e,Ao)' 



(6) 



where /„(e, 0) is the first principles value for the T — 
K contribution to the inelastic current of tunneling holes 
with energy e, and Aq = 0.6 eV, £5,7 — Sp ^ —6 
eV(including band bending) are the parameters at T = 

We determine the T dependent width. At, by fitting 
R calculated using In{e,T) and 7t to the measured val- 
ues at Vb = — 7 V. The fitted values of At are plotted as 
function of T in Fig. 4, and solid curves in Fig. 2 show 
the corresponding T dependent values of R as function 
of / for Vb = -5 V and Vb = -7 V. At Vb = -5 V 
the dashed curves and the solid curves coincide, and the 
T dependence of R mainly stems from 7t, whilst there 
is a strong additional T dependence of R from At at 
Vb = —7 V. This behaviour can be understood using 
R ^ {12/lT)^'^ and taking into account the T depen- 
dence of I2 from Eq. ra). Using Eq. (H) it is observed 
that I2 oc At for |Vb| <C \ep — escrl and I2 oc Aj^^ for 
|Vb| > |e|. - £5a|. At Vb = -7 V this gives a T de- 
pendence of i?(300)/i?(380) - (1.3 X 1.135)6-5 ~ 10^ 
and i?(300)/i?(450) ^ 10^ in good agreement with the 
measured data. At Vb = — 5 there is a weak T depen- 
dence of R upon At since V^l < \£f~ ^5<y\ causes I2 to 
be nearly independent of T and the T dependence of R 
therefore mainly stems from 7t. The solid curves plot- 
ted in Fig. 3 as obtained from the present model are in 
good agreement with the experimental data for the whole 
bias range from —10 V to —4 V. It is observed that the 
curves become more flat with increasing temperatures 
in good agreement with experimental data. At positive 
bias inelastic scattering with the Si-H 6cr* is the dom- 
inant desorption mechanism for |Vb| < 3. In this bias 
regime |Vd| <S l^f;' — SQa* | ~ 6 eV and at this condition 
the main T dependence of R must come from 7t in good 
agreement with Ref. g]. 

In order to investigate the origin of the T-dependence 
of At we model At by p2] 



At = An 



a(e 



e/kT 



1)- 



(7) 



From a least-squares fit to the data in Fig. 4 we obtain 
Ao = 0.49 eV, a = 0.36 eV, and 9 = 402 K. The T 



dependence of At is quite strong and almost linear for 
T > 300 K with the slope a/Q - 1 meV/K. The mea- 
sured dependence upon T of Aph for silicon is in the 
range 0.1- 0.4 meV/K [g2|,g3|, and a major part of the 
T dependence must therefore come from Ac. This de- 
pendence must be related to a change in the electronic 
coupling between the adsorbate and the surface upon ex- 
citations of phonons with frequency TiVi ~ Q. However, 
the fit in Eq. (4) is quite insensitive to the value of 8 and 
from the present data we can only state that TiQ, < 0.1 
eV. It is therefore not possible to distinguish whether the 
T dependence of A is from Si phonons or Si-H transverse 
vibrations. 

In summary we have presented experimental data of 
the desorption rate of STM induced desorption of H from 
the Si(100)-H(2xl) surface as a function of temperature, 
and found a decreasing desorption rate as function of 
the temperature. The desorption mechanism is related 
to vibrational heating of the H atom by holes inelastic 
scattering with the Si-H ba hole resonance, and the tem- 
perature dependence is explained by both a decreasing 
vibrational life time of the Si-H stretch mode and a de- 
creasing electronic life time of the Si-H 5(t hole resonance 
with increasing temperature. 
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FIG. 1. Scanning tunneling micrograph of exposed silicon 
dangling bonds as a result of hydrogen (H) desorption on a H 
terminated Si(100)-(2xl) reconstructed surface of the n-type 
sample at a substrate temperature of 530 K and a line scan 
with a speed of 2 nm/s, a sample bias (Vb) of —7 V and a tun- 
nel current (/) of 8.0 nA. The filled-state image was recorded 
at Vb = -1.7 V and I = 0.2 nA. 
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FIG. 2. STM induced desorption rate (R) as function of 
tunnel current (/) for temperatures 300 K, 380 K and 450 K 
for sample biases of (a) —7 V and (b) —5 V. Curves show theo- 
retical calculations of R using only a temperature dependent 
vibrational lifetime (dashed) and combined with a tempera- 
ture dependence of the electronic lifetime(solid). 
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FIG. 3. The tunnel current (7dcs) as function of sam- 
ple bias (Vb) which gives rise to 50% desorption of H from 
a Si(100)-H(2xl) surface along a line scan of 2 nm/s by an 
STM tip. Results are plotted for temperatures between 300 
K, 380 K, 450 K, 570 K, and 610 K. Curves show theoretical 
calculations(see caption of Fig. 2). 
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FIG. 4. Squares show temperature dependent values of A 
obtained by fitting theoretical desorption rates to the exper- 
imental values. The solid curve shows a least-squares fit of 
Eq. (m) to the data points. 



